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' The theory of tidal interactions between an accretion disc 
, and embedded companion has been elaborated in a series 
of papers by Lin & Papaloizou (1979a & b, 1986a & b) (see 
also works by Hourigan & Ward 1984, Goldreich & Tremaine 
1980). These papers were primarily concerned with the early 
solar nebula, and the question of orbital migration of the 
proto planets, in particular of the giant planets. An an- 
nular gap may be formed in the accretion disc due to the 
action of the tidal forces from the planet, and the formation 
of such a gap leads to a coupling between the orbit of the 
planet and the viscous evolution of the disc. In the existing 
work the satellites were light compared to the local disc, so 
they migrated on a local viscous timescale. There was little 
modification of the disc properties outside the gap. In this 
case, therefore, the presence of the satellite modifies the disc 
spectral energy distribution only through the removal of a 
narrow annulus of emitting material around the satellite. 

In this paper we explore the effects of a satellite that is 
massive compared with the local disc. Table 1 gives details 
of two astrophysical situations in which this could be the 
case. The inertia of the satellite prevents it from migrating 
on the local viscous timescale. As a result, disc material is 
banked up upstream. Downstream of the satellite, material 
is able to accrete onto the central object, so that the mass 
of the disc interior to the satellite decreases with time. The 
satellite therefore causes global changes in the disc structure, 
and enhances the relative contribution of material upstream 
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to the spectral energy distribution emitted by the disc. 

We consider the coupled evolution of a satellite-disc sys- 
tem in the limit that the disc extends only upstream of the 
satellite. Depletion of the inner disc means that the influ- 
ence of the inner disc on the satellite dynamics is relatively 
minor. The neglect of the inner disc implies either that the 
disc formed from material of higher specific angular momen- 
tum than the satellite (as in the case of a collapse scenario 
where the satellite formed from low angular momentum ma- 
terial, with the disc being formed subsequently from higher 
angular momentum), or else that the satellite is sufficiently 
massive for the inner disc to have largely drained onto the 
central object. In any case, as the satellite migrates radi- 
ally, it becomes relatively more massive than the local disc, 
so we anticipate that our calculation of the evolution of the 
satellite and disc structure is insensitive to the omission of 
the inner disc. The resultant spectral energy distribution at 
high frequencies may be more sensitive to this assumption. 

The structure of the paper is as follows. In Section 2 we 
review the process of gap formation and in Section 3 we de- 
scribe an analytic model for the evolution of the disc and the 
orbit of the satellite. In Section 4 we describe two numeri- 
cal approximations to the same problem, and compare their 
results with the analytic predictions. In the final section 
we discuss some astrophysical applications of our models, in 
active galactic nuclei, and in protostellar discs. In the AGN 
case we show that presence of massive compact objects in 
the disc would modulate the emission on a timescale of up 
to 10 5 years. In the case of a protostellar disc, we consider a 
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satellite which is a massive giant planet or sub-dwarf. The 
long timescales for the inward transit (~ 10 5 years) implies 
that such satellites would be expected to cause a substan- 
tial fraction of T Tauri stars to show spectra markedly differ- 
ent from those generated by standard steady-state accretion 
discs. 



2 THE EVOLUTION OF COUPLED 
DISC-SATELLITE SYSTEMS 

Suppose the satellite has a mass M s and its orbit is at a 
radius r s . Consider a 'collision' between the satellite and a 
parcel of gas with mass fi. The effect of the tidal torque of 
the planet on the disc may be thought of by analogy with 
dynamical friction: the fluid at r > r s moves slower than the 
satellite, and hence exerts a drag on it, reducing its angular 
momentum and tending to drive its orbit to smaller radius. 
The fluid at r < r s does the opposite. The effect on the disc 
of the back reaction to the friction is to push the fluid away 
from the planet; to decrease the specific angular momentum 
of the fluid at small radii, and to increase it at large radii. 
The response of the disc is determined by a competition 
between these tidal perturbations from the satellite, tending 
to open up an annular gap around the satellite, and pressure 
and viscous forces in the disc, tending to close the gap. 

We will derive an expression for the torque using the 
analogy with dynamical friction. This is essentially the 
calculation done by Lin & Papaloizou (1979a), but in the 
present case, since our purpose is illustrative, we omit con- 
stants of order unity. Lin & Papaloizou (1986a) examine the 
validity of the impulse approximation numerically. (See also 
Lin & Papaloizou 1986b for a discussion of the same point.) 
Consider a single collision with a fluid particle at an impact 
parameter p in the impulse approximation (e.g. Binney & 
Tremaine 1987). If the relative velocity of the satellite and 
fluid is v, then the angular momentum transferred to the 
satellite is 



Ai~fi 



G 2 M 2 r s 
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where G is the gravitational constant. Averaging over such 
impacts, we find that the rate of change of L is 



dL 
dt 

which we write as 



AL — vdp 
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dL 
dt 



EATr s dp 



with the torque per unit mass of the disc given by 

G 2 M 2 
At = — r-. 



(2) 
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Writing q for the mass ratio of the satellite to central object, 
M s /M„, and expanding v in a Taylor series about p = we 
have 



A T = g 2 fi 2 r s 2 ( — 



(5) 



where fi s is the angular velocity of the satellite, assuming 
a Keplerian rotation law and circular orbits for the satellite 
and disc. 

The specific viscous torque near the edge of the gap, 
assuming that the radial scale length ~ A, is 



fir 

A v ~u — 



(6) 



where v is the viscosity in the disc and fi is the local angular 
velocity. The radial extension of the gap, A, can be evalu- 
ated approximately by equating A v and At with p ~ A: 



A 



fir 2 2 

— q • 



(7) 



Note that in the limit v — > equation (7) correctly predicts 
that A — > oo, since in this case there is no equilibrium gap: 
the satellite continually eats into the disc, at a slower rate 
as the gap grows. Following Lin & Papaloizou, we define a 
dimensionless measure of the gap width 



A 



(8) 



which is the ratio of the viscous to orbital timescales, mul- 
tiplied by the square of the mass ratio. 

Two criteria need to be fulfilled in order for the satellite 
to succeed in evacuating an annular gap in the disc. Firstly 
there is the requirement that the gap be not re-filled by 
radial pressure gradients, which means that the disc vertical 
scale height, H , must be less than the width of the gap 



A > H. 



(9) 



Secondly, there is the requirement that the effect of the satel- 
lite is to transfer angular momentum, rather than to accrete 
from the disc. This means that the gap width should be 
greater than the Roche radius of the satellite: 



A > R L = q 1/3 r. 
This reduces to 



q > 



fir 2 



(10) 



(11) 



i.e. q must exceed the local inverse Reynolds number. In 
the relatively cool and inviscid discs which satisfy these cri- 
teria, the satellite sits in an annular gap in the disc of width 
A, exchanging angular momentum with the disc material 
through the action of tidal torques concentrated near the 
gap's edge. 

In principle it is possible to arrange interior and exterior 
disc material such that the net torque acting on the satellite 
is zero, but such a situation will not persist once the disc has 
undergone significant viscous evolution. The depletion of 
material interior to the satellite by accretion will result in an 
imbalance between the torques experienced by the satellite. 
The net angular momentum transfer is outward from the 
satellite, which therefore spirals in. 

Let us define the term 'reference disc' to mean the 
steady state disc in the absence of the satellite. We will 
refer to quantities in the reference disc with the subscript 
0. Now we define a dimensionless parameter measuring the 
mass of the satellite compared to the reference disc 



B 



47rS r s 2 



(12) 
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where Eo is the surface density in the reference disc. A 
relatively light satellite, in the sense that B ^ 1, does not 
substantially modify the surface density of the disc outside 
the gap. It behaves like a representative fluid element in the 
disc, and therefore its orbit evolves radially on the viscous 
timescale (Lin and Papaloizou 1986b). 

In this paper we are interested in the opposite case, 
where B ^ 1. The disc is unable to push the satellite radially 
inwards on its own viscous timescale and material therefore 
forms a reservoir upstream of the gap. Let us define 



(13) 



where v is the viscosity in the disc. Suppose that the thermal 
and hydrostatic equilibrium relation between v and E has 
the power law form given implicitly by 



cab 

Locj r , 



(14) 



where a and b are constants. The largest contrast in E, over 
the reference value Eo, allowed by condition (11), is then 



Ea 

E^ 



Or" 



2 \ a/(l-a) 



(15) 



where vq is the corresponding reference value of v. If the 
surface density behind the satellite rises to E s ~ E a then 
the gap will be overwhelmed by viscous diffusion. 

In a Shakura-Sunyaev (1973) disc, where the viscosity 
is given by 



(16) 



with a constant, the criteria (9) and (10) reduce to, respec- 
tively, 



< 



H_ 

r 

and 
H 



< 



1/2 



1/5 



(17) 



(18) 



If (17) is violated, the gap will be closed by pressure; if 
(18) is violated, the gap will be closed by accretion onto 
the satellite. In most cases of astrophysical interest H/r 
is an increasing function of r in the reference disc, so once 
the satellite opens up a gap it persists, as long as the disc 
structure is not modified by the satellite (i.e. if B > 1). 
If the disc structure changes as a result of the presence of 
the satellite, and H/r becomes a decreasing function of r, 
the values of q and a determine whether it is accretion or 
pressure which eventually closes the gap. If 



1/2 



1/5 



(19) 



then if the gap closes, accretion will be responsible. This 
is illustrated in Figure 1 in which we have defined three 
important radii: r is the radius at which the criterion (10) is 
satisfied; tb is the radius at which 5 = 1; and r c is the radius 
at which (10) is no longer satisfied in the perturbed disc. 
Figure 1 is entirely schematic, and no particular viscosity 
law has been used to calculate H/r. When r s < r a gap is 
opened in the reference disc, and when r s < tb the disc is 
banked up behind the satellite. 




Figure 1. The aspect ratio H/r of the disc as a function of r. 
The dot-dashed line is the reference disc. The long-dashed line 
shows the effect of a satellite in banking up the disc. The lower 
and upper short-dashed lines are the critical values for closure of 
the gap by accretion and pressure respectively. The solid line is 
a possible trajectory of H/r immediately behind the satellite as 
it migrates radially. 

When r s < r c or r s > r the satellite cannot maintain an 
annular gap in the disc, so that the analysis involved in the 
present work (which assumes tidal coupling between satel- 
lite and disc across such a gap) is invalid. Knowledge of the 
evolution of a satellite embedded in a disc flow requires two- 
dimensional hydrodynamic calculations beyond the scope of 
the present work. Thus the following remarks are necessar- 
ily speculative. Once the satellite reaches the radius r c , we 
would expect that the inward flow of material as the gap is 
overwhelmed would result in a burst of accretion luminosity. 
The fraction of this material that is accreted by the satellite, 
rather than by the central object, is however not clear. It is 
also not clear whether the satellite would remain embedded 
in the disc thereafter or whether the overflow of a small frac- 
tion of the dammed up mass would allow the gap to re-open. 
In the latter case the satellite might continue to spiral in, 
maintaining a state of marginal stability against gap closure, 
with intermittent overflow of the material upstream. 

The maximum energy which could be released when the 
gap is closed is 



E b 



Gtf«S c 



(20) 



This corresponds to accretion of the mass AM ~ E c r;? by 
the central object of mass M* and radius r*. 

The expectations that a satellite would move as a rep- 
resentative fluid element in the case B > 1 was confirmed 
by Lin and Papaloizou (1986b) using one-dimensional nu- 
merical calculations. In their calculation, the usual viscous 
diffusion equation for the disc was supplemented by an addi- 
tional tidal torque term (equation 5). The main motivation 
for these calculations was to evaluate migration timescales 
for the proto- Jupiter. The question of substantial damming 
up of the disc by a satellite considerably more massive than 
the local disc was not extensively investigated. In the fol- 
lowing section we consider this question in some detail. 
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Table 1. The approximate values of relevant physical param- 
eters in a selection of astrophysical situations. 

M* M s r B t B 

AGN 10 6 M@ 10 3 M@ ~ 10 4 GM/c 2 10 5 y 

T Tauri star 1M & 10~ 2 M & ~ 100_R @ 10 5 y 

Af* is the central mass; M s is the satellite mass; r B is the radius 
within which the satellite is more massive than the local disc; 
i_B is the radial migration timescale of the satellite at r B . 



3 AN ANALYTIC MODEL FOR 
SATELLITE-DISC EVOLUTION 

In this section we first describe the characteristics of the ref- 
erence disc. We then develop a model for the surface density 
of the disc when influenced and enhanced by a satellite with 
B < 1 (equation 12). The slow radial drift of the satel- 
lite enables viscous torques in the disc to communicate the 
presence of the satellite to larger radii. We characterise the 
upstream structure of the disc by dividing it into two sec- 
tions: an outer region where the influence of the satellite has 
not yet been felt, and one in the vicinity of the satellite. The 
boundary between these regions is denoted rt- The outer re- 
gion, rt < r, is assumed to be in a quasi-steady state with 
a mass flux determined by the fuelling rate at infinity. The 
inner region, r s < r < rt, is assumed to be in a quasi-steady 
state the properties of which are determined by the angular 
momentum transfer between the satellite and the inner edge 
of the disc. This results in a lower mass flux in the inner 
region , on account of the damming effect of the satellite. 
The limit of rt — > oo is analytic and explicit. 

3.1 The reference disc 

In a thin, Keplerian accretion disc with no external torques, 
/ is determined by 



9E _ 3d_ L dr'f 

dt r dr dr 

(e.g. Pringle 1981), and so in a steady state 



d L dr 2 f 
— r 2 i- 

dr dr 



0. 



(21) 



(22) 



In a steady-state disc we assume an outer boundary 
condition 



Ldr'f Mo 
lr2 —> — = T— = Too, r -> oo, 
or o7T 



(23) 



where Mo is the mass flux through the disc. At the inner 
boundary, r», we assume the torque on the inner edge of the 
disc vanishes, whence 



/ = 0, r = r*. 

The general solution to (22) is 
f = X + ^ 

f 2 



(24) 



(25) 



where X and Y are constants. With the boundary condi- 
tions (23) and (24) we find that 



f — fo — foo — fa 



(26) 



In the limit r* — > 0, which is adequate for our purposes, 

fo = foo- 




100. 



Figure 2. The surface density S s behind the satellite as a func- 
tion of r s (short dashed line), in the analytic model with rt oo. 
The constants of the model are chosen so that B = 0.04 at r = 1. 
The solid lines are snapshots of the surface density S. Also shown 
is the reference disc (long dashed line). 



3.2 The satellite 

We consider the situation in which the disc interior to r s is 
empty (E = 0). Such a situation would arise if the disc was 
formed in this way, but would also be approximately realised 
in the limit of a sufficiently massive satellite: the satellite 
acts as a dam, allowing the inner disc to drain away in a time 
much shorter than the migration timescale of the planet. 

We assume that the tidal coupling between the satellite 
and the disc is confined to a very narrow region at the in- 
ner edge of the gap. In this case equation (21) adequately 
describes the evolution of /. The outer boundary condi- 
tion remains the same, on the assumption that far enough 
upstream of the satellite the disc must approach the undis- 
turbed solution. The inner boundary condition is deter- 
mined by supposing that a) the satellite and disc edge move 
with a common velocity, and b) there is overall conservation 
of angular momentum. The radial velocity in the disc is 
given by conservation of mass as 



3 dr'f 

Er2 <9r 



(27) 



Evaluating this at r = r s , and equating it to the velocity 
of the satellite, v s , provides the inner boundary condition. 
(Provided A remains considerably smaller than r s , then v 
and v s cannot substantially differ.) 
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Angular momentum conservation is guaranteed by 
equation (21) everywhere except at the inner edge of the 
gap. The viscous torque on the edge of the disc is 

C !-//<> (28) 

assuming a Keplerian rotation law. The rate of change of 
the angular momentum of the satellite is 

1 

2" 



^-M s r s 2 Q s 
dt 



-M s r 2 Q s — 



(29) 



again assuming a Keplerian rotation law. Conservation of 
angular momentum is thus expressed as 



!-/ / <> + -.1/ / <> — = 0, 

2 r s 

so the value of v in the boundary condition (27) is 
67rr s / s 

" s = mT' 

We write equation (31) as 



-B 



fo to 



phere 



to 



2r_i_ 
3 i/o 



(30) 



(31) 



(32) 



(33) 



is the viscous timescale in the undisturbed disc. 

In reality the disc will never reach a global steady state, 
since v s is always non-zero. But, for small B, the satellite 
moves slowly compared to the local viscous timescale, so we 
make one more approximation: that the disc is able to reach 
a steady state for all values of r s . In the outer region of the 
disc, rt < r, we suppose that it has not had time to adjust 
to the presence of the satellite, and thus that 



f = fo=U 



rt < r. 



(34) 



In the inner region, r s < r < rt, we solve (22) with the 
boundary conditions (34) at r = rt, and a combination of 
(27) and (32) at r = r s . Using (14) and (25), we find that 
the required solution to (22) is 



(35) 



/ = /a + (/» " /a) - 



where 37r/ a < Mo is the mass flux through the inner disc. 
Equating (34) and (35) at r = rt we find that 



/a + (/, " /a) ( - ) 2 = U 



with 



-l/(l + a) 



(36) 



(37) 



Once rt is specified, equations (36) and (37) are sufficient to 
determine f s and / a , given r s . 

The model assumes that the disc is in a steady state 
everywhere apart from at rt . At rt the mass fluxes in the 
inner and outer disc are unbalanced, so mass is not conserved 
unless rt moves at just the right rate. The evolution of rt 
would be specified by equations (36) and (37), together with 
an integro-differential equation that ensured conservation of 
mass. We do not develop the latter in this paper, but retain 
the notion of an inner and outer disc for the purpose of 
qualitatively understanding our numerical simulations. 



3.3 The limit r t -> oo 

In the limit rt — > oo, the model becomes explicit and an- 
alytic. We expect that this limit should be approached as 
B — > 0, in which the satellite moves very slowly, and is able 
to influence the surface density very far upstream. In this 
limit we find that 



/» = foB- 



-l/(l + a) 



(38) 



This is an approximate description of the 'banking up' of 
the surface density behind the satellite. The surface density 
just behind the satellite is given by 

E 5~ a/(1 + a) . 



E s 



(39) 



Provided that a > we are reassured that the smaller the 
value of B, the greater is the surface density E s . The surface 
density can thus be enhanced so that E s > Eo. This is 
illustrated in Figure 2, in which we have used a viscosity 
law (equation 14) with a = 0.65 and b = —0.65, which is 
appropriate to quiescent protostellar discs (Clarke, Lin & 
Pringle 1990). In all the quantitative work in this paper, 
and in all the figures (except Figure 1), we use the same 
viscosity law. 

Combining (38) and (32) we can find v s explicitly in 
terms of r s : 



-B 



a/(l + a) Jj_ 

to 



(40) 



Integrating this equation we find the evolution of r s with 
time 



n \l- C jB« /(1+a ^ 1/C 



(41) 



where ti and 5; are the values of to and B evaluated at 
r = r;, and 



2 + b 
1 + a ' 



(42) 



Comparing equation (40) to the radial velocity at r s in the 
undisturbed disc 



vo 



r s 
'To' 

we find that 

v s = voB /v ' . 



(43) 
(44) 



For B < 1, a > 0, b > —2, we see that v s < vq. Our 
assumption of a quasi-steady state must break down when 
v s <^ vo — but we expect v s to be at most vo, when B ~ 1. 



4 NUMERICAL MODELS 

4.1 Infinitessimal boundary layer 

To test the steady-state assumption in the analytic model, 
we now drop that assumption, but retain the inner bound- 
ary condition (equations 27 and 31), and then solve equation 
(21) numerically. The righthand side of equation (21) is dif- 
ferenced in a standard implicit scheme, and a Lax average 
is employed for the lefthand side (e.g. Press etal. 1992). 
An iterative procedure is employed for solving the implicit 
equations, as described for example by Potter (1987). A uni- 
form grid in the variable r 1 / 2 is employed, with 6400 points 
in the range r = r; n = 0.25 to r = r ou t = 1024. The initial 
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Figure 3. The result of dropping the steady-state assumption: 
snapshots of f(y) with B = 0.04 at r = 1 (solid curves). The 
short-dashed curve shows the analytic prediction for the right- 
most snapshot with rt = 30; the long-dashed curve shows the 
analytic prediction with rt oo. 



satellite. A more massive satellite causes a larger surface 
density enhancement and is able to influence the disc out to 
a larger radius. The analytic model is in each case able to 
fit the numerical results adequately, once rt is specified. 

4.2 Explicit torque 

As a second numerical comparison with the analytic model 
we drop the assumption of an infinitessimal boundary layer, 
and include the tidal torque from the satellite explicitly. 
Thus equation (5) appears as a source term in the diffusion 
equation (21): 



df Id, _ 
-J- = 3r 2 

dt r dr \ d 



dr~f A T E 



(45) 



(Lin & Papaloizou 1986b). For practical purposes, At has 
to be softened on a scale which is smaller than A. We use 
the prescription of Lin & Papaloizou (1986b), but since E 
vanishes in the gap anyway, the details of the softening are 
not important once the gap has formed. 

Angular momentum conservation is ensured by setting 
the rate of change of the angular momentum of the planet 
equal to the integrated contributions of the tidal torque ex- 
erted by the disc: 
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Figure 4. As Figure 3 but with B = 0.01 at r = 1. The 
short-dashed curve in this case is with rt = 200 for the right-most 
snapshot; the long-dashed curve shows the analytic prediction 
with rt oo. 

radius of the satellite orbit is r = 6, and the inner boundary 
condition is applied at the nearest gridpoint upstream. At 
every timestep, the new location of the satellite is computed 
from equation (31), and if it crosses a gridpoint, the inner 
boundary condition is applied at that gridpoint. 

In Figure 3 we show two snapshots of the disc structure 
for the case B = 0.04 at r = 1; we plot f j fo as a func- 
tion of y = r -1 / 2 . The figure illustrates that the disc is 
indeed well represented by a sequence of quasi-equilibrium 
states — f(y) can be approximated by two straight line sec- 
tions. The radius of the intersection of these straight lines, 
rt, does not change much as the satellite migrates radially. 
(We have checked the independence of rt on the position 
of the outer boundary by carrying out additional calcula- 
tions with r ou t = 256, in which the results were the same.) 
Figure 4 (B = 0.01 at r = 1) shows the case of a heavier 



The numerical resolution has to be high for small A 
(small gap width). Here we used 12800 grid elements in 
the radial range r £ [0.25,256]. The results of the explicit 
torque calculation are shown in Figure 5. In common with 
Figure 3, we see that / is approximately described by two 
straight line segments as in the analytic model, and that the 
location of rt does not change much with time. This illus- 
trates that, as expected, the form of the torque distribution 
has negligible effect on the disc structure at large distances 
upstream. The disc structure close to the satellite is some- 
what different when the distributed torque is included. The 
disc is effectively obeying the boundary condition (equation 
31) at a radius where the tidal torque is maximum. But the 
tidal torque does not immediately drop to zero outside this 
radius, and the surface density can continue to rise. As a 
result, the peak surface density is greater than that implied 
by equation (31). The distributed torque model predicts 
higher peak densities than the infinitessimal boundary layer 
model. 

In Figure 6 we compare the results of using two different 
values of A (equation 8). A smaller value (Figure 6a) leads 
to a narrower gap, and a sharper turn over of the surface 
density profile at the gap edge, but does not significantly 
affect the value of the surface density at the turn over. We 
conclude that the gap width (A ~ A 1 ^ 3 r s ) has a minor ef- 
fect on the surface density enhancement, which is controlled 
largely by the ratio of the disc mass to the mass of the satel- 
lite (B). 



4.3 Summary 

Figure 7 shows a comparison of the evolution of r s for the 
numerical models with the analytic model. In the case of 
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Figure 5. The result of including an explicit tidal torque: the 
explicit torque model with B = 0.04 at r = 1 and A = 0.25. The 
dashed curve shows the analytic prediction for the right-most 
snapshot with rt oo. 

the infinitessimal boundary layer model, the surface density 
enhancement is smaller than predicted by rt — > oo. The 
surface density enhancement determines the radial evolution 
of the satellite through equation (31), so the evolution of 
r s is slower than the analytic prediction. In the case of the 
explicit torque model it is slower because the radius at which 
equation (31) is satisfied is closer to the satellite than the 
peak surface density. 

The numerical results show that the disc has a quasi- 
equilibrium structure with rt varying slowly. The limit rt — > 
oo, in which all the equations are analytic and explicit should 
be approached in the limit of a massive satellite. Since the 
numerical calculations are considerably more expensive in 
effort and computational resources, the advantages of the 
analytic approach are obvious. 

5 DISCUSSION 

We have derived analytic expressions for the surface density 
enhancement in a disc which is modified by the presence of 



Figure 7. The evolution of r s as a function of time (in units 
of the viscous time to at r = 1), for B = 0.04 at r = 1. The 
long-dashed line is the explicit torque model; the short-dashed 
line is the infinitessimal boundary layer model; and the solid line 
is the analytic model with rt oo. 

a massive satellite at its inner edge (equations 12 and 39). 
We have verified the accuracy of these expressions, to order 
unity, using two levels of numerical approximation. 

A satellite initially at large radius, where B > 1, moves 
likes a representative fluid element and scarcely modifies the 
structure of the disc. Provided B decreases with decreasing 
radius (b > —2, equation 14), such a satellite may eventually 
reach the radius, rs, at which 5 = 1. The disc surface den- 
sity is subsequently enhanced behind the satellite because it 
is relatively massive compared with the local disc (B < 1). 

In Figure 8 we compare the spectral energy distribu- 
tion, F v , in a disc dammed up behind a satellite with that 
generated by a reference disc with a zero torque boundary 
condition at its inner edge. In all cases we assume that the 
disc radiates everywhere as a black body with local effec- 
tive temperature fixed by the local viscous dissipation rate. 
When the satellite is at radii considerably larger than the 
inner edge of the reference disc, the assumed absence of ma- 
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Figure 8. The evolution of the spectral energy density F v from 
the accretion disc, as a satellite drifts radially inwards. The 
dashed line shows F v for the reference disc, with inner radius 
r* = 0.5, r ou t = 200 and a = —b = 3/5. The solid lines show 
F v from the banked-up disc, which is assumed to be empty inside 
the radius of the orbit of the satellite. From left to right, the 
solid curves are for r s = 20, 10 and 5. In Each case rt oo and 
B = 0.01 at r = 1. 

terial interior to the satellite implies that the spectral energy 
distribution is deficient in higher frequency radiation com- 
pared with the reference disc. (This effect would be more 
marked if the contribution from boundary layer emission 
were included — in the case that the reference disc accreted 
onto a stellar surface at its inner edge). As the satellite pro- 
gresses to smaller radii, the emission is enhanced over a wide 
range of frequencies: the energy lost by the inward spiralling 
satellite is radiated through the disc, and the enhanced emis- 
sion from just upstream of the satellite compensates for the 
absence of material interior to it. The spectrum of a classic 
steady state accretion disc may be characterised as a power 
law F v oc v 1 / 3 over a wide dynamic range. When a satel- 
lite with B < 1 is present this changes to F v oc j/ 5 / 7 (see 
Appendix for derivation). 

We consider two applications of astrophysical interest: 
a disc around a pre-main sequence (T Tauri) star with a 
sub-dwarf satellite; and a disc in an Active Galactic Nucleus 
where the satellite is either of stellar mass or is a compact 
object of mass ~ 10 3 Mq (see Syer, Clarke and Rees 1991 
and Lin 1994 for arguments relating to the possible presence 
of such objects in AGN discs). The T Tauri model uses a 
central mass of IMq, viscous alpha parameter 10 -3 and an 
accretion rate of 10 -r MQy -1 (details of the disc structure 
in Clarke, Lin and Pringle 1990). The AGN model uses a 
central mass of 10 6 Mq, viscous alpha parameter 10 -1 and 
an accretion rate of one tenth the Eddington rate (details of 
the disc structure in Clarke 1987, 1988). 

For a stellar mass object in an AGN disc, the disc over- 
whelms the gap created by the satellite at effectively all 
radii. When the mass of the satellite is 10 3 Mq, and in the 
T Tauri model the gap opening criteria are satisfied and it is 
meaningful to calculate tb, as listed in Table 1. A massive 
compact satellite in an AGN disc would bank up the disc 
behind it for r s < re ~ 10 4 times the radius of the central 
object (at this radius the disc is cool and unionised). In the 



pre-main sequence case, a sub-dwarf of mass 1O" 2 M would 
bank up the surface density in the disc for r s < tb ~ IOORq. 
The time taken for the companion to spiral in from tb is of 
order the disc's viscous timescale at tb', t s , the time to spiral 
in from a smaller radius is greater than the local reference 
viscous timescale, to, because the inertia of the satellite is 
large compared with the local disc. At smaller radii the 
ratio t s /to increases. However, Ib provides a measure of 
the interval of time over which the satellite is able to signifi- 
cantly modify the structure and spectral energy distribution 
of the disc. Thus a massive black hole in an AGN disc would 
give rise to deviations from classic steady-state disc emission 
during its transit of the inner disc. It may also affect the 
continuum emission, on the assumption that the continuum 
arises from physical processes near the central mass, and 
therefore relies on the disc to transport fuel to small radii. 
A realistic calculation of the spectral energy distribution 
would also have to include the remnant disc interior to the 
satellite, and would require knowledge of the physical pro- 
cesses that give rise to the continuum emission. We find that 
for a large range of AGN disc parameters, the disc is unable 
to close the annular gap created by the companion and that 
the banking up of the disc behind the satellite (and associ- 
ated increase in luminosity) continues until the companion 
reaches the last stable orbit and is swallowed by the central 
black hole. At this point the spectral energy distribution of 
the disc reverts to that of the usual steady state disc on a 
timescale comparable with Ib- 

In the pre-main sequence case the corresponding value 
of ts is long (~ 10 5 years) so that the disc structure and 
spectrum would deviate from its steady state values over a 
considerable period. Although, as above, the calculation of 
the spectral energy distribution should also contain a com- 
ponent from the disc inside the satellite (a contribution that 
becomes increasingly unimportant once the satellite is well 
inside tb), Figure 8 gives some indication of the types of 
spectral energy distributions that might be generated. Thus, 
contrary to previous attempts to model the effect of satel- 
lites on disc emission merely in terms of "missing emission" 
from annular gaps (e.g. Marsh and Mahoney 1992), we find 
that satellites massive compared with the local disc produce 
global modifications in the spectral energy distribution: the 
disc emission is enhanced upstream of the satellite and this 
produces a variety of signatures in the resultant spectral en- 
ergy distribution. While the satellite is at relatively large 
radius, the increased strength of the outer disc produces a 
hump in emission at low frequencies. This is in line with 
the observed tendency of T Tauri stars to exhibit too much 
power at low frequencies compared with standard steady- 
state disc models. Once the satellite is very close to the star, 
its effect is to steepen the spectral energy distribution from 
the classic i/ 1 / 3 to a i/ 5 / 7 power law. This effect would be 
masked observationally by the contribution from the bound- 
ary layer between the depleted disc and the central object. 

In the case of a satellite in a T Tauri disc we find that 
the disc is never able to overwhelm the gap created by the 
satellite, provided it remains cool and largely unionised; if 
this is the case, the satellite continues to spiral inwards until 
it is either enveloped or tidally disrupted by the central star. 
Another possibility is that the temperature in the disc up- 
stream of the satellite reaches the threshold for the hydrogen 
ionisation; if this is the case then the associated rapid rise in 



Satellites in Discs 9 



the viscosity could enable the disc to overwhelm the gap and 
flow rapidly inwards past the satellite. Such a scenario may 
offer a mechanism for the generation of FU Orionis outbursts 
in pre-main sequence stars (Clarke & Syer in preparation). 
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where k is Boltzmann's constant. Thus the spectrum of the 
disc as a whole is 



font 

F v cx / B V (T) 2-kt dr 



(49) 



where r; n and r ou t are the inner and outer radii of the disc, 
respectively. 

Suppose that 

T = n(^fY (50) 

in which case for hv <C kT* we can rewrite equation (49) as 
Iout z 1+2 'f>dz 



,2//3 3-2/(3 
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(51) 



where x = hvjkT and x oll t = hv/kT oll f For frequencies 
such that kT oll t <C hp <C kT* we may take x oll t 3> 1 and 
hence F v cx j/ 3_2//3 . 

In the reference disc fi = 3/4 and we recover the classic 
result that F v oc v 1 / 3 . When the disc is truncated by a 
satellite we take r; n = r s , and assume that the disc interior 
to the satellite is empty. For B <C 1, equations (35) and 
(38) imply that 



-i/(i + a) / r. 



1/2 
\ /" / 

Thus we find that f3 = 7/8 and 



(52) 



(53) 



Thus, the spectrum of the disc is steeper than the reference 
disc, and the normalisation increases as B decreases (i.e. as 
r s decreases) for a > — 1. 



This paper has been produced using the Blackwell Scientific 
Publications Tp^X macros. 
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THE EMITTED SPECTRUM OF 



The response of the disc to the tidal forces from the satellite 
causes addtional viscous dissipation in the disc. If the inner 
disc is evacuated, then there is also a deficit of high fre- 
quency radiation compared with the reference disc. Here we 
calculate the spectral energy density of the radiation from 
the disc in our analytic model. We show that the former ef- 
fect dominates at late times, and the latter dominates when 
the satellite first begins to stem the radial flow in the disc. 

On the assumption that the disc is optically thick, at 
each radius there is an effective temperature of the emitted 
radiatio T given by 

/ 



crl cx 



(47) 



where a is the Stefan-Boltzmann constant (e.g Pringle 
1981). The spectrum from each annular element of the disc 
is the blackbody spectrum 



B V (T) cx 



exp(hv / kT) 



(48) 



